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O ■ ABSTRACT 

The disk around AB Aur was imaged and resolved at 24.6 /xm using the Cooled 
\0 , Mid-Infrared Camera and Spectrometer on the 8.2m Subaru Telescope. The 

gaussian full- width at half-maximum of the source size is estimated to be 90 ± 6 

AU, indicating that the disk extends further out at 24.6 /im than at shorter wave- 
. lengths. In order to interpret the extended 24.6 /im image, we consider a disk with 

Q_i. a reduced surface density within a boundary radius Re, which is motivated by 

Q ! radio observations that suggest a reduced inner region within about 100 AU from 

^— > ! the star. Introducing the surface density reduction factor fc for the inner disk, 

c^ , we determine that the best match with the observed radial intensity profile at 

24.6 /im is achieved with Rc=S8 AU and /c=0.01. We suggest that the extended 

> 

f^ I ^ Based on data collected at Subaru Telescope, which is operated by the National Astronomical Observa- 

OO ■ tory of Japan. 
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emission at 24.6 yum is due to the enhanced emission from a wall-hke structure 
at the boundary radius (the inner edge of the outer disk), which is caused by a 
jump in the surface density at R^. Such reduced inner disk and geometrically 
thick outer disk structure can also explain the more point-like nature at shorter 
wavelengths. We also note that this disk geometry is quahtatively similar to a 
pre-transitional disk, suggesting that the AB Aur disk is in a pre-transitional 
disk phase. 



Subject headings: protoplanetary disks 
sequence 



circumstellar matter — stars: pre-main 



Introduction 



The evolution of protoplanetary disks has attracted attention because it is presumably 
a site of the on-going planet formation. Mostly by analyzing spectral energy distributions 
(SEDs), many tran sitional disks whose iii ner region is devoid of small dust grains have been 
detected. Recently, lEspaillat et al.l ( 120071 ) introduced a new sub-class called pre-transitional 
disks with definite near-infrared (NIR) excesses that indicate the presence of an inner disk 
separated from an optically thick outer disk. They suggested that such a disk is forming 
a gap and when the inner disk dissipates, it will become a transitional disk with an inner 
cavity. In both cases, the presence of a wall-like structure on the inner edge of the outer disk 
can explain the excess mid-infrared (MIR) emission seen in the SED. For interr aediate-mass 



young stars such as Her big Ae/Be stars, so me objects such as HD 135344B ( IGrady et al. 



20091 ) and HD 100546 (JBenisty et al.l l2010l ) are thought to be in the pre-transitional or 
transitional disk phase. 



AB Aur (d=144;^^7pc; AO Ve, Ivan den Ancker et al.l (119971 )) has been regarded as a pro- 
totypical Herbig Ae/Be star and is o ne of the best-studie d Herbig Ae stars, with a mass of 
2.4±0.2 Mq and an age of 4±1 Myr (JDeWarf et al.l 120031 ). Previous observations at various 
wavelengths revealed complex circumstellar structures at various scales. Optical and NIR 
coronagraphic direct images showed an extended circumstellar nebulosity of the envelop e 
and spiral structures in the disk in scattered light (JGrady et al.lll999t iFukagawa et al.ll2004j ). 
Recent coronagra phic polarimetric imaging observations revealed a dust annulus at a ra- 
dius of ~100AU JOppenheimer et al.l l2008riPerrin etaPbood : iHashimoto et al.lboiOJ ) with 
a possible fluctuation in the azimuthal direction. Millimeter and sub-mm interferometric 
observations sho wed that the gas and dust have be en depleted in the central region ( r <100 
AU) of the disk Jpietu et al.lboosi : llin et al.ll2006l ). 
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At MIR wavelength s, many direc t imag ing: and interferometric studies have been per- 
formed toward AB Aur. I Chen fc Jural (120031 ) imaged AB Aur at 11.7 and 18.7 yum using the 
Keck I telescope, and resolved it at 18.7 /im. They found that models with a disk surrounded 
by an envelope fit their data better than flared disk models. Applying a deconvolution tech- 
nique to the 2 0.5 Atm data obtained at the 3.6 m Canada- France- Hawaii- Telescope (CF HT) , 
Pantin et al.l (120051 ) discovered a ring-like structure at 280 AU. iMariiias et al.l (120061 ) pre- 
sented 11.6 and 18.5 /im images using the Gemini North telescope; however, they could not 
find the ringlike structure reported by iPantin et al.l ( l2005l ) . They derived the intrinsic source 
full-width half-maximum (FWHM) to be 17±4AUat 11.6 /im and 22±5 AU at 18.5 /im, as- 
suming a gaussian brightness distribution. iLiu et al.l (l2005l . l2007l ) also resolved the disk at 
10.3 /xm by nulling interferometric observations and derived the diameter of the disk to be 
24-30 AU. 



In this paper, we present results of the first imaging observations of AB Aur at 24.6 /im 
using the 8.2m Subaru telescope. At 24.6 /im, the point-spread function (PSF) is stable 
compared to that at shorter wavelengths because of the larger Fried length, which enables us 
to observe small extended structures with high reliability. In addition, it allows us to trace 
the cooler part of the disk. By constructing a simple model, we discuss the disk structure of 
AB Aur as well as implications for the nature of the AB Aur disk. 



2. Observations and Data Reduction 



We m ade imaging observa t ions with the Cooled Mid-Infrared Cam era and Spectrometer 
(COMICS; iKataza et aPboOol : lokamoto et"aDl2003[ ISako et al.l[2003l l on the 8.2 m Subaru 
Telescope on Mauna Kea, Hawaii on 2003 October 12. AB Aur was observed using the 
Q24.5-OLD filter (A^„,_=24.6/im, AA=1.9/im). The plate-scale of COMICS is 0.13 arcsec 
per pixel. The chopping throw was 10 arcsec and the position angle (PA) of the chopping 
direction was 30 degree. The chopping frequency was 0.45 Hz, and the total integration 
time was 2005 s. Just before and after AB Aur observations, we observed a Tau as the PSF 
reference and photometric standard star. The total integration times were 401 s and 251 s, 
respectively. 

For data reduction, we employed a shift-and-add method to improve the spatial reso- 
lution. The imaging data consist of 0.983 s on-source integration frames. First, the ther- 
mal background and the dark current signals were removed through the subtraction of the 
chopped pair frames. The object was sufficiently bright to be recognized even in 0.983 s in- 
tegration chop-subtracted frames, so we searched for the peak position by fitting a Gaussian 
with the source profile. Then we shifted the frames (in order to align the peak position) and 
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summed them, but we excluded the frames whose gaussian FWHM deviated more than 1 a 
from the mean value. This rejection of lower quality data resulted in an effective integration 
time of 1333 s for AB Aur. The same procedure was applied to the data of a Tau taken 
before and after AB Aur observations, and the effective integration times were reduced to 
327 s and 209 s, resp ectively. Flux calibr ation was performed using the a Tau template 



spectrum provided by lEngelke et al.l (120061 ). Standard aperture photometry was performed 



and the resultant flux density of AB Aur was 43.0 ± 3.4 Jy. 

3. Results 

3.1. Source size of AB Aur at 24.6 /im 

The 24.6 /xm images of AB Aur and a Tau are shown in FiglH It is clear that AB 
Aur is extended beyond the PSF at this wavelength. The disk emission appears symmetric 
and we could not find clear evidence of the radial profile depending on the azimuthal angle 
within the uncertainty. This is probably beca use the disk around AB Aur is nearly face- 



on (inclination ~ 30°; iFukagawa et al.l l2004j ). The FWHMs of shift-and-added images 



of AB Aur, as well as a Tau before and after observations of AB Aur were 0.82" ±0.03", 
0.647" ±0.007" and 0.644" ±0.004", respectively. This error is the standard deviation of the 
selected frames. The FWHM of the PSF reference [a Tau) is comparable to the predicted 
value of the diffraction-limited performance of the telescope at 24.6 /im. We estimate a rough 
source size of AB Aur by comparing the observed AB Aur image with a circular gaussian 
image convolved with the observed PSF image. The best-match gaussian FWHM is 90±6 
AU. This value is larger than that at shorter wavelengths, since the gaussian FWHMs of the 



AB A ur disk are repoted to be 17±4 AU at 11.6 /im and 22±5 AU at 18.5 /im ( iMariiias et al. 



20061). 



3.2. Comparison with a simple model 

Gaussian brightness distributions are used for the sake of simplicity and they are not a 
good representation of the true intrinsic radial brightness profile of the disk. Thus, it would 
be useful to compare the observed AB Aur image with a disk model for the investigation 
of physical properties of the disk. We can calculate the model brightness distribution of 
a disk if the temperature and density structure of the disk is given. H ere we derive the 



disk structure based on a 1±1D radiative transfer model similar to that in lDuUemond et al. 



( I2OO2I ): the stellar radiation transfer is solved by the so-called grazing angle approximation, 
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but the radiative transfer of the dust thermal radiation is solved exactly along the disk 
vertical axis with a variable Eddington factor method. The dust and gas are assumed to 
be thermally and dynamically coupled, and the vertical density structure is determined by 
the hydrostatic equilibrium that is maintained consistent with the temperature structure. 
Scattering is not taken into account. When the disk is optically thick, the model predicts 
a flared geometry and reproduces the spectrum of protoplanetary disks. When the disk 
becomes optically thin, the grazing angle approximation does not work. In such cases, we 
treat the stellar radiation transfer by a 2D ray-tracing method. This modification is found to 
work in the present calcula tion. We assume dust consisting of spherical astronomical silicate 
grains (JDraine fc Led 119841) and the dust size distribut i on giv en by n{a) oc a~^'^, where a is 
the grain radius ( iMathis et al.l 119771 ). As iPerrin et al.l (|2009[ ) discussed that the maximum 
grain size is about l/xm, we set the minimum and maximum grain sizes to be 0.005 /im 
and 2.5 /xm, respectively. The stellar parameters used are: the e ffective temperature T ^, 
= 9,500 K, the radius R, = 2.5 R©, and the mass M, = 2.4 M© Jvan Boekel et al.ll2005h . 
We assume a power-law surface density S = Si(r/AU)~^ between ri^ and rout, where p=l 



(JHartmann et al.lll998l ). Si = 70 g cm 



0.5 AU and r^ 



1000 AU. The surface 



density is the total (gas + dust) one and we also assume a gas-to-dust mass ratio of 100. We 
chose El and rout so as to reproduce the observed far- infrared to sub-mm/mm flux densities 
and rin is the radius where the dust temperature becomes about 1,400 K. In addition, we 
introduce the boundary radius Re that separates the inner and outer regions of the disk 
because the inner region (within 70 — 110 AU of the c entral star) is reported to show a 



deficiency in the mm/sub-mm dust continuum emission (jPietu et al.l 12005 



Lin et al.l 120061 ) 



which is attributed to a decrease in the surface density of dust grains smaller than mm- 
size in the inner region. To account for this information in our disk modeling, we consider 
that the surface density is reduced by a scaling factor fc of 0.1, 0.01, and 0.001 for r < Re 
compared to the power-law surface density of the outer disk. The explored parameter space 
is summarized in Table 1. This reduction of the inner surface density makes the inner 
region much thinner than the outer disk, and as a result, a wall- like structure appears at the 
boundary radius. The upper part of the wall is directly illumina ted by the central star; thu s 
its temperature is enhanced like a surface super- heated layer (IChiang fc GoldreichI 119971 ). 
On the other hand, the lower pa rt of the wall remains cool because it is shadowed by the 
inner region (JMulders et al.ll2010l ). In practice, we calculate the disk structures of the inner 
and outer regions separately, and connect them at the boundary radius Re. The temperature 
of the illuminated upper part of the wall (i.e. inner edge of the outer region) is set to be 
the same as the surface layer at the wall. In addition, the same treatment is adopted at the 
most inner edge of the disk which is also directly illumi nated by the centra l star. We need 
2D radiative transfer calculations to be self-consistent (JMulders et al.l 120101 ). which will be 
done in the future. 
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On the basis of the derived disk structure, we calculate the model brightness distribution 
for various wavelengths and spectrum by a ray-tra cing method along the observer's line of 
sight. The inclination angle is assumed to be 30° (JFukagawa et al.ll2004f ). To compare the 
model brightness distribution with our observation, we create a model image where the 
model brightness distribution is convolved with the observed PSF image shown in Fig. 1. 
Finally, we find that the best match with the observed radial profile is obtained with the 
scaling factor of the inner disk /^ = 0.01 and the boundary radius Re = 88 AU. This is 
comparable to the inner radius derived by iPietu et al.l ( l2005l ) (~110 AU for dust and ~70 
AU for CO). Fig. 2 shows the best-fit model brightness distribution at 25 /im, the model 
image for Subaru/COMICS and its peak- normalized radial profile. Structural parameters 
of the best-fit model such as the radial distributions of the surface density and height, are 
shown in Fig. 3. The surface height was determined by the height where the stellar Planck 
mean opacity became unity towards the central star. Th e spectrum of the best- f it model also 
matches with that of A B Aur from 10 /xm to ~ 1 mm (JHelou fc Walkerlll988t iMeeus et al. 



200ll : lAcke et al.ll2004i and references therein) as shown in Fig. 4. For wavelengths shorter 
than 10 /xm, our model spectrum underestimates the obs erved spectrum of AB A ur. This 
may be due to the emiss ion from the inne r-most gas disk ( jTannirkulam et al.ll2008l ) or from 
the puffed-up inner rim (jNatta et al.ll200ll ). neither of which is considered in our model. The 
contribution of these inner structures to the 25 fim flux density is negligible; the spectrum at 
1-8 yum can be well fitted by an additional blackbody of temperature 1,400 K (the upper solid 
line) and this blackbody component emits 1.4 Jy at 25 /im, only 3% of the 25 /im Infrared 
Astronomical Satellite (IRAS) flux density of 48 Jy (see Fig. 4). 



4. Discussion 
4.1. Comparison with other observations and validity of the vi^all-like structure 



A reduction of the surface density inside the boundary radius {Re ~ 88 AU) and the 
resultant wall-like structure at Re can explain the extended emission we observed at 24.6 /xm. 
Additionally, this interpretation is consistent with the mm/sub-mm interferometric observa- 
tions of AB Aur. Furthermore, it can explain the properties of the AB Aur disk that appear 
more point-like at shorter MIR wavelengths. In order to determine the extent of the AB Aur 
disk at these shorter MIR wavelengths from our best-fit model, we calculated 11.4 and 18 
/im model brightness distributions of this model, and convolved them with a gaussian model 
PSF comparable to the diffraction limit of an 8.2m telescope. Then we measured FWHMs 
of the PSF convolved model images and derived gaussian FWHMs of 8 AU at 11.4 /zm and 
21 AU at 18 /im by using the quadratic subtraction method. The value at 18 /xm is con- 
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sistent with previous observational resul ts by iMariiias et al.l (120061 ) and that at 11.4 /xm is 
slightly narrower than the observations ( IMariiias et al.ll2006l ). but the general trend is fine. 
The more compact nature of the AB Aur disk at the shorter MIR wavelengths relative to 
25 /im is due to not only the radial temperature distribution but also the temperature of 
the wall-like structure (96 K). This temperature is too low for dust to radiate efficiently at 
shorter wavelengths, resulting in little change of the appearance of the disk at shorter wave- 
lengths. A similar situation is also seen in another Herbig Ae/Be sta r HP 142527, which is 
only marginally extended at 10 /im, but clearly extended at 24.5 /im (JFujiwara et al.ll2006l ). 



Another piece of collateral evidence for the wall-like structure at ~ 88 AU is the agree- 
ment betwe en the temperature (93K) of the blackbody component derived from the spectrum 
of AB Aur (IMeeus et al.ll2001) and the temperature (96K) of the wall at 88 AU in our best- 
fit disk model. 



Meeus et al. 



( 120011 ) suggested that this temperature might refiect the dust 
temperature at the onset of disk fiaring. Their suggestion is qualitatively similar to ours in 
the sense that the wall in our model can be interpreted as the onset of the outer disk fiaring. 

Furthermore, the annulus of dust at ~100 AU is observed i n the NIR polarimetri c 



coronagraphic images (jOppenheimer et al.l l2008l : iPerrin et al.l l2009l : JHashimoto et al.l |2010| ) . 
Inside the d ust annulus of the NIR scattered hght, the gap is reported at a radius of 0.5" ~0.6" 
(72~86AU; JHashimoto et al.ll2010l ). This also supports the onset of the disk fiaring and the 
presence of a wall-like structure at ~88 AU in the AB Aur disk. 



4.2. The AB Aur disk: a pre-transitional disk? 

Combined with the relatively young age of the central star and the pristine environment 
as indicated by the global envelope structure, AB A ur has been though t to be a prototypical 



example of newly formed Herbig Ae/Be star (e.g. JGrady et al.lll999l ). However, the inner 
region of the AB Aur di sk possesses little gas and dust a s revealed by mm/sub- mm inter- 
ferometric observations (jPietu et al.ll2005l : JLin et al.ll2006l ). indicating the evolved nature of 
the disk. Such an a pparently paradoxical situation may have led to the difficulty faced by 
Chen fc Jural (120031 ) in explaining the AB Aur properties with a simple fiared disk model. 
A reduction in the inner disk surface density such as that used in our model can solve this 
problem. This situation is qualitati vely similar to pre-transitional disks introduced for some 
T Tauri stars ( jEspaillat et al.ll2007| ) in the sense that an inner disk separated from the opti- 
cally thick outer disk is present. We propose that the AB Aur disk is also a pre-transitional 
disk. 



Various mechanisms for the clearing of the inner region have been suggested such as 



Alexander et al 



dust coagulation (e.g. iDullemond fc Dominikl 120051). inside-out evacuatio n by the magnetro- 



tatio nal instability (e . g. IChiang fc Murray- Clayl 1200 71 : ISuzuki et al 



20061 ). dynamical effects of companions (e.g. 



20091) . photoeyaporat ion 



Varniere et al.ll2006al Jbl). 



(e.g. 

and so on. In the case of AB Aur, it is not clear which mechanism is the most plausible; 
however, the prese nce of planetary companions has been suggested, which may have cause d 
the inner hole (e.g. iPietu et al.ll2005t iJang-Condell fc Kuchneiil2010l : iHashimoto et al.ll2010l ). 
The pre-transitional disk nature of the AB Aur disk may reflect ongoing planet formation. 
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Fig. 1. — Shift-and-added image of AB Aur (top) and a Tau as a PSF reference (bottom) at 
24.6 /im. The surface brightness is given in Jy/arcsec^ and scaled from to the peak value 
of each image. North is up and east is to the left for the AB Aur image; however, the PSF 
image is aligned to match the same instrument rotation angle as the AB Aur image for a 
similar optical configuration. 
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Table 1. Range of Parameters Explored in Model Grid 



E scaling factor /c 



Boundary radius Re [AU] 



0.1 
0.01 

0.001 



75, 81, 88, 95, 102, 110, 119, 128 

65, 70, 75, 81, 88, 95, 102, 110, 119 

33, 38, 44, 52, 60, 65, 70, 75, 81, 88, 95, 102 



model brightness 
distribution 




radial profile 

n AB Aur (obs.) 
°n A PSF (obs.) 
model*PSF 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 
radius [arcsec] 



Fig. 2. — Plots comparing the observations and the best-fit model (boundary radius Re = 88 
AU and scaling factor fc = 0.01). The left panel represents the model brightness distribution 
at 25 /im in Jy/arcsec^. The middle panel shows the model image where the model bright- 
ness distribution is binned to the COMICS plate-scale (0.13" /pixel) and convolved with the 
observed PSF shown in Fig. l; additionally, the im age was rotated so that the disk major- 
axis position- angle was 58° (JFukagawa et al.ll2004j ). The right panel shows peak- normalized 
azimuthally-averaged radial profile plots of the observed AB Aur image (squares), the PSF 
image (triangles), and the binned model image convolved with the PSF shown in the middle 
panel (diamonds); horizontal axis represents radius in arcseconds. 
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Fig. 3. — (a) Radial distribution of the summed gas and dust surface density for the best-fit 
model (boundary radius Re = 88 AU and scaling factor /c = 0.01). (b) Radial distribution 
of the surface height Hs normalized by R for the best-fit model. The surface height was 
determined by the height where the stellar Planck mean opacity became unity along the line 
of sight to the central star. A wall-like structure is present at the boundary radius Re- 
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Fig. 4. — Spectrum of AB Aur observed and tha t of our best - fit mo del. Observed data are 
from lHelou fc Walked ( 19881 ): iMeeus et al.l ( 200ll ): lAcke et al.l ( 2004J) and references therein. 
Redde ning correctioii i s app lied using Av=0.5 mag ( JAcke et al.ll2004j ) assuming the extinction 
law of ICardelli et al.l ( 119891 ). The stellar emission is shown in dotted line. The dot-dashed 
line represents emission from the dust disk of our best-fit model. The dashed line is for the 
additional 1400K blackbody needed to fit the observed spectrum, which is not considered in 
our best-fit model. The lower solid line is a sum of the emission from the star and the dust 
disk, while the upper solid line is a sum of the star, disk and 1400K blackbody emission. 



